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The behavior of solitons in integrable theories is strongly con- 
strained by the integrability of the theory, that is by the existence 



of an infinite number of conserved quantities that these theories are 
known to possess. As a result the soliton scattering of such theories are 
expected to be trivial (with no change of direction, velocity or shape). 
In this paper we present an extended review on soliton scattering of 
two spatial dimensional integrable systems which have been derived 
as dimensional reductions of the self-dual Yang-Mills-Higgs equations 
and whose scattering properties are highly non-trivial. 
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g : 1 Introduction 

CO 

Solitons which scatter in a non-trivial way and which occur as solutions 
of planar integrable systems are going to be presented. Initially solitons 
were introduced by mathematicians to describe lumps of energy stable to 
\ perturbations which do not change either their velocity or their shape when 

colliding with each other. They have been observed experimentally both as 
waves on shallow water and in laser pulses in fibre-optic cables, among other 
places. Solitons also play an important role in various models of subatomic 
particles (quantum field theories). In recent literature all sorts of localized 
energy configurations have been called solitons which travel without assuming 
stability of the shape or velocity or a simple behavior in collision. 

In one space dimension integrable systems occur when dispersion effects 
are exactly balance by nonlinearities. However, in more than one dimensions 
the definition of an integrable model is related with the existence of infinite 
number of conserved quantities; of Painleve property; of two linear equations 
(so-called the Lax pair) which compatibility condition give the non-linear 
soliton equation; of inverse-scattering transform; even of multi-soliton solu- 
tions (see, for example, [I]). 
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An interesting problem is the investigation of the scattering properties 
of two or more solitons colliding. In some known models with non-trivial 
topology the collision of two solitons is inelastic (some radiation is emitted) 
and non-trivial, ie a head-on collision results in 90° scattering (see, for exam- 
ple, |2] and references therein). On the other hand, the solitons of integrable 
models interact trivially in the sense that they pass through each other with 
no lasting change in velocity, shape or direction. Some examples in 2+1 
dimensions are the Kadomtsev-Petviashvili the Konopelchenko- Rogers 
[3], the Davey-Stewrtson |3] and the integrable chiral equation c 6j. Until 
recently, non-trivial scattering of solitons occurs mostly in non-integrable 
models which is far from simple. The issue discussed here is whether such 
type of scattering can occur in integrable models too. There are some limited 
examples of integrable models where soliton dynamics is non-trivial. In 1+1 
dimensions there are many models that non-trivial soliton-like solutions (cf. 
Ref. [7]) exists like the boomeron solutions jH] which are solitons with time- 
dependent velocities. In 2+1 dimensions there are the dromion solutions 
of the Davey-Stewartson equation which decay exponentially in both spatial 
coordinates and interact in a non-trivial matter and the soliton solu- 
tions of the Kadomtsev-Petciashvilli equations whose scattering properties 
are highly non-trivial jTTj. 

Recently in families of soliton solutions have been constructed 

for planar integrable systems (using analytical methods) with different types 
of scattering behavior. This happens since the solitons of the systems have 
internal degrees of freedom that although they determine their space orien- 
tation, they do not change the energy density and they are important in 
understanding the system evolution. Therefore, solitons can interact either 
trivially or non-trivially depending on the orientation of these internal pa- 
rameters and on the values of the impact parameters defined as the distance 
of closest of approach between their centers in the absence of interaction. 

In what follows we will review the approach based on the Riemann-Hiblert 
problem with zeros, first presented in |I2]-P3], i n order to construct soliton 
solutions with non-trivial scattering for two spatial dimensional integrable 
models defined in anti-de Sitter and Minkowski space, respectively. Both 
models are equivalent to the self-dual Yang-Mills-Higgs equations reduced 
from 2+2 to 2+1 dimensions while the soliton solutions are related to hyper- 
bolic and Euclidean Bogomolny-Prasad-Sommerfield (BPS) monopoles. 
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2 Self-dual Yang- Mills- Higgs Equations 



Static BPS monopoles are solutions of non-linear elliptic partial differential 
equations on some three-dimensional Riemannian manifold. Most work on 
monopoles has dealt with the case when this manifold is the Euclidean space 
R 3 since then the equations are integrable and sophisticated geometrical 
techniques like twistor theory can be applied. Note that the introduction of 
time-dependence destroys the integrability. In addition, the monopole equa- 
tions on hyperbolic space H 3 are also integrable and often hyperbolic 
monopoles turn out to be easier to study than the Euclidean ones as first 
studied by Atiyah and explicitly shown in jTHj. In fact as it has been rig- 
orously established in |T7j, in the limit as the curvature of hyperbolic space 
tends to zero the Euclidean monopoles are recovered. 

Let us first consider an integrable system JH] which is related to hyper- 
bolic monopoles and which is obtained from replacing the positive definite 
hyperbolic space H 3 by a Lorentzian version, the so-called anti-de Sitter 
space. Thus, the self-dual Yang-Mills-Higgs equations are of the form 

D^ = ^= gij e jkl F kl (1) 
2^ 

where the Yang-Mills-Higgs fields take values on a three-dimensional Rie- 
mannian manifold {M.) with gauge group SU(2). In particular, (for 
k — 0, 1, 2) is the su(2)-valued gauge potential, F^ = diAj — djAi + [Aj, Aj] 
is the field strength and $ = ^(x^) is the sw(2)-valued Higgs field; while 
) represent the local coordinates on JA. The action of the 
covariant derivative Di = di + Aj on $ is: D$> = 8$ + [Ai, $]. Equation ((TJ) 
for constant curvature is integrable in the sense that a Lax pair exists. 

Note that the solutions of which can be described in terms of holomor- 
phic vector bundles or in terms of rational functions correspond to Euclidean 
or hyperbolic BPS monopoles when (Ai,g) is the Euclidean R 3 or the hyper- 
bolic H 3 space, respectively. Also, both of the models presented here can be 
obtained from (Q) as dimensional reductions of the four- dimensional self-dual 
Yang-Mills-Higgs equations for appropriate gauge choices [T]. 
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3 The Anti-de Sitter Model 



Currently a great deal of attention has been focused on anti-de Sitter spaces 
since they arise naturally in black holes and p-branes. For the case of Yang- 
Mills theory with M = 4 supersymmetries and a large number of colors it has 
been conjectured that gauge strings are the same as the fundamental strings 
but moving in a particular curved space: the product of five-dimensional anti- 
de Sitter space and a five sphere JH]- Then, using Poincare coordinates the 
anti-de Sitter solutions play the role of classical sources for the boundary field 
correlators, as shown in .20, ; while extensions of the corresponding statements 
can be applied to gravity theories, like black holes which arise in anti-de Sitter 
backgrounds. 

By definition the (2+l)-dimensional anti-de Sitter space is the universal 
covering space of the hyperboloid Tt defined by the equation 

U 2 + V 2 - X 2 - Y 2 = 1 (2) 

with metric 

ds 2 = -dU 2 - dV 2 + dX 2 + dY 2 . (3) 
By parametrizing the hyperboloid 7i by 

U = sec p cos 9 

V = sec p sin 9 
X = tan p cos <fi 

Y = tan p sin <fi (4) 
for p G [0, 7r/2), the corresponding metric simplifies to 

ds 2 = sec 2 p [-d9 2 + dp 2 + sin 2 p rf0 2 ) . (5) 

The spacetime contains closed timelike curves due to the periodicity of 9 |2*T] . 
In fact (p, 0) correspond to polar coordinates and 9 G R being the time; while 
anti-de Sitter space (as a manifold) is the product of an open spatial disc 
with 9 and curvature equal to minus six. Null spacelike infinity X consists of 
the timelike cylinder p = n/2 and this surface is never reached by timelike 
geodesies. 
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Figure 1: The Penrose diagram of anti-de Sitter space. The boundary of 
anti-de Sitter is the boundary of the cylinder. 



If the Poincare coordinates (r, x, t) for r > are defined as 

1 



x 



t 



U + X 
Y 

U + X 
-V 



U + X 

the metric simplifies to the following form 



ds 2 



r' 2 (-dt 2 + dr 2 + dx 2 



(6) 



(7) 



Note that the Poincare coordinates cover a small part of anti-de Sitter space 
which correspond to half of the hyperboloid 7i for U + X > as shown in 
Figure ^ The surface r = is part of infinity I. 
Consider the set of linear equations 



rD r - 2 (A - u)D u - $] * 
A — u 



A — u 
2D V + D 




0. 



(8) 



Here A £ (C and (r, u, v) are the Poincare coordinates (where u = x + 1 and 
v = x — t) while the gauge fields ($, A r , A u , A v ) are 2x2 trace-free matrices 
depending only on (r,u,v) and ^(A, r, u, t) is a unimodular 2x2 matrix 
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function satisfying the reality condition \E'(A)\[ / (A)^ = / (where t denotes 
the complex conjugate transpose). The system (jSJ) is overdetermined and in 
order for a solution \I/ to exist the following integrability conditions need to 
be satisfied 

D u <$> = rF ur 

D r $ = -2rF uv (9) 

which correspond to the self-dual Yang-Mills-Higgs equations (Q) defined in 
the 2+1 anti-de Sitter space. The gauge and Higgs fields in terms of the 
function \I/ can be obtained from the Lax pair (jHJ) due to the boundary 
conditions. Note that, as A — > oo the function \1/ goes to the identity matrix 
and the system (jHJ) implies that 

A u = 0, A r = - $. (10) 
r 

On the other hand, for A = and using (JTUJl the rest of the gauge fields are 
defined as 

T 

$ = --J r J —uJ u J 

Av = j r j-i-j v j-i (ii) 
zr 

where J(r,u,v) = \I/(A = 0,r,u,v). As a result the first equation of the 
system © is automatically satisfied (due to the specific gauge choice). 

Recently, Ward has shown that holomorphic vector bundles V over Q 
determine multi-soliton solutions of in anti-de Sitter space via the usual 
Penrose transform. This way a five-parameter family of soliton solutions 
can be obtained, in a similar way as for flat spacetime |22]- Later, more 
solutions of equations (0) were obtained by Zhou (221 121] using Darboux 
transformations with constant and variable spectral parameter. In what 
follows, we use the Riemann problem with zeros to construct families of 
soliton solutions and observe the occurrence of different types of scattering 
behavior. More precisely, we present families of multi-soliton solutions with 
trivial and non-trivial scattering. 
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3.1 Baby Monopoles 

In this section, we illustrate the construction of time-dependent solutions 
related to hyperbolic monopoles. In particular, families of soliton solutions 
of the self-dual Yang-Mills-Higgs equations defined in the (2+l)-dimensional 
anti-de Sitter space are constructed and their dynamics is studied in some 
detail. 

The integrable nature of means that there is a variety of methods for 
constructing solutions. Here, we indicate a general method for constructing 
soliton solutions of Q which is a variant of that presented in Ref . [22] • Using 
the standard method of Riemann problem with zeros, in order to construct 
the multi-soliton solutions, we assume that the function \l/ of the system (jSJ) 
has the simple form in A 

* = ' + £t^ (12) 

where are 2x2 matrices independent of A and n corresponds to the soliton 
number. The components of the matrix Mk are given in terms of a rational 
function fkiyJk) = dk^k + Ck of the complex variable: = v — r 2 (/!*; — u)~ x . 
Here a^, and fik are complex constants which determine the size, position 
and velocity of the A;-th solitons. Remark The rational dependence of the 
solutions \1/ follows (directly) when the inverse spectral theory is considered. 
In fact in j2S] (for flat spacetime), it was shown from the Cauchy problem 
that the spectral data is a function of a parameter similar to u^- 

By considering the unitarity condition of the eigenfunction \1/ and keeping 
in mind that the gauge fields are independent of the spectral parameter A it 
can be shown that the matrix Mk has the form 

n 

Mk = Y,(^ 1 ) M m l a m k b (13) 
i=i 

with r _1 the inverse of 

a=l 

and holomorphic functions of Uk of the form m k = (m^m^) = (1, 
The Yang-Mills-Higgs fields ($, A r , A v , A u ) can then be read off from (fTUl - 
ITTj) and they automatically satisfy while the corresponding solitons are 
spatially localized since $ — > at spatial infinity. 
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Figure 2: A two-soliton configuration with trivial scattering at different 
times. 
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By way of example, let us look at the special case where /ii = 1 + i, 
fi 2 = 2z, ai = 2, a 2 = 1, C\ = 5 and c 2 = —10. Figure 121 represents snap- 
shots of the positive definite gauge quantity (— tr$ 2 ) at different times. The 
corresponding solution consists of two solitons which travel towards infinity 
(r = 0) and bounce back while their sizes change as they move. 

Note that the Riemann problem with zeros method was first introduced 
by Zakharov and his collaborators i n his pioneer work of applying spec- 
tral theory to generate soliton solutions of integrable equations. Later this 
approach has been applied to obtain the monopole solutions of the four- 
dimensional self-dual Yang-Mills-Higgs equations [27]. However, as it can 
been observed from (|12|) . the Riemann problem with zeros method assumes 
that the parameters fik are distinct and also Jx^ ^ fii for all (k, I). In what 
follows examples are given of two generalizations of these constructions: one 
involving higher-order poles in /!& and the other where flk = fii- 



<0 High-order poles in 



Firstly, let us look at an example in which the function \1/ has a double 
pole in A and no others. In this case, \I/ has the form 

- ,+ £iAf (16) 

where Rk are 2x2 matrices independent of A. Then \1/ corresponds to a 
solution of (jHJ) if and only if it factorizes as ^1] 

^V-T^JDl^jV-Tx^ltr) (16) 

for some two vectors q and p. One way to derive the structure of these 
vectors is to take the formula (|12|) for n = 2, set fi\ = /i + e, \ii = \i — e 
and fxiuji) = f{u x ) + sh(u x ), / 2 (w 2 ) = /(w 2 ) - eh(u 2 ), with / and h being 
rational function of one variable. In the limit e — > the two vectors q and p 
can be obtained and are of the form: 

q = (l + | / | 2 )(lJ) + (^-^)^^^ + /^(/,-l), 

P = (1,/) (17) 
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where / and h are rational functions of u = v — r 2 (fi — u)^ 1 while /' denotes 
the derivative of / with respect to its argument. In this case, the solution 
depends on the parameter fi and on the two arbitrary functions / and h. 
Note that the constraint faiuv) — — » as e — > has to be imposed in 
order for the resulting solution ^ to be smooth for all (r, u, v) which is the 
case here. 

In order to illustrate the above family of solutions two simple cases are 
going to be examined by giving specific values to the parameters /i, f(u) and 
h(u). 

• Let us study the simple case where \i = i, f{uj) = uj and h{uj) = 0. 
Then, the quantity — tr$ 2 simplifies to 



- tr$ 2 - 32^ ^ + x2 ~ t2 + 1 ^ + 4t ') {^1 + ^ 2 -t 2 -l) 2 + 4x 2 
f((r 2 + x 2 - t 2 ) 2 + 1 + 2t 2 + 2x 2 ) 2 + 4r 4 l 2 



;is) 



which is time reversible. The corresponding time- dependent solution 
is a travelling ring-like soliton configuration which for negative t, goes 
towards spatial infinity (r = 0); approaches it at t — and then bounces 
back at positive t while the soliton size deforms, as shown in Figure El 
Ring structures occur in the soliton scattering of many non-integrable 
planar systems and are approximations of multi-solitons 



• Next, we investigate the solution which corresponds to a non-trivial 
soliton scattering as presented in Figure 0] The solution has been 
obtained by (jl6H17|) for the following values of the parameters: fi = i, 
f{uj) = uj and h(u) = uj 4 while the picture consists of two different-sized 
solitons. For large (negative) t, the gauge quantity — tr$ 2 is peaked at 
two points, forms a lump at t — and then two solitons emerge, for 
large (positive) t, with the small one been shifted to the left. 

These solutions can also been obtained using Uhlenbeck's construction 
In this approach the function \1/ is assumed to be a product of factors 
which when substituted in the Lax pair (jSJ) the problem simplifies to finding 
solutions of simple first-order partial differential equations (201 • Both meth- 
ods can be extended to derive solutions where the function \1/ has higher order 
pole in A (and no others). Then, \1/ can be written as a product of three (or 
more) factors with three (or more) arbitrary vectors (for more details, see 
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;ure 3: A stationary soliton configuration at different times. 
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Figure 4: A non-trivial dynamics of two solitons in anti-de Sitter space. 



B3|). 

^> Case where pk = 

Secondly, let us construct a large family of solutions which correspond 
to the case where fik = One way of proceeding is to take the solution 
(I12j) with n — 2, put fii = fi + e, \ii = ft — e and take the limit e — > 0. In 
order for the resulting \l/ to be smooth it is necessary to take fi(ui) = f(u)i), 
12(^2) = — 1//(^2) — £h(u 2 ) where / and h are rational functions of one 
variable. On taking the limit we obtain a solution \I/ of the form 

n 1 <S> rn 1 n 2 <8> m 2 . . 

V = I+- + — — 19 

A — fi A — fi 

where n k , m k for k — 1,2 are complex valued two vector functions of the 
form 

.1 fl f\ ™2 



m 1 = (l,/), m 2 = (-/,!) 



2(^-/2) /2(1 + |/| 2 ) -(H)w\/m lf 

nV 4(l + |/| 2 ) 2 -(//-//) 2 | W | 2 V(/i-p)^ -2(1 + |/| 2 )/U 2t 



(20) 



with 



So we generate a solution which depends on the parameter \i and the two 
arbitrary rational functions / = f(ui) and h = h(u). 

In Figure 03 we represent snapshots of the solution (JT§j) for fj, — i, f — 
uj and h = Q. The configuration consists of two solitons with non-trivial 
scattering behavior. Again, the quantity — tr$ 2 is peaked at two points, for 
(negative) t, which are still distinct at t — and then two shifted (compared 
to the initial ones at t — —2.5) solitons emerge, for (positive) t. Throughout 
the time-evolution their sizes change. Note that, the scattering solutions 
belong to a large family since / and h can be taken to be any rational 
functions of uj. 

Remark: The extension of the obtained classical solutions in the whole 
anti-de Sitter space, ie using the coordinates (p,0,4>), is unambiguous. For 
example, the simplest solution which corresponds to the one soliton (first 
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Figure 5: A non-trivial dynamics of two solitons in anti-de Sitter space. 



derived in [TH]) given by (|12|) for n = 1, //i = i and /i = uj\ implies that 

8r 4 



tr$ 2 



[(r 2 + - t 2 ) 2 + 2x 2 + 2t 2 + l] 2 
2 cos 4 p 

2 n _ 9\1 



(cos 2 p — 2 



(22) 



which means that the positive definite quantity — tr $ 2 is independent of the 
variables (8, d>). 



4 The Minkowski Model 

The integrable SU{2) chiral model introduced by Ward |H] is described by 
the field equation 

d v (j-'duj) - d x (J-^J) = (23) 

where the chiral field takes values in: J 6 SU(2) while u = \{t + y) and 
v = ^(t — y). Note that the integrable chiral equation is related to the Jarvis 
equation studied in [T£|. Once more the equation can be obtained from the 
self-dual Yang-Mills-Higgs equation given by (0) for specific gauge choice 
(further details in [3*T]). Then for specific boundary conditions spherical 
symmetric monopole solutions could be derived using the harmonic map 
ansatz as studied in [16J. 

Here the choice of the boundary conditions follows from the chiral-equation 
form rather than the gauge-theory form, ie 

J=Jo + ;Ji(*)+o(^), r^oo (24) 

for x + iy = re td . Here Jo is a constant matrix and Ji depends only on 9 (no 
time dependence). 

The model possess a conserved energy density given by jHj 

£ = -^tr jf + (J^oU) 2 + (J- I d y j) 2 ^j (25) 

which is a positive-defined functional of the chiral field. 

The integrable equation (}2*3*|) can be written as the compatibility condition 
of the system 

(Xd x -d u )iP = At 

(Xd v -d x )iP = B^j (26) 
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where A G (D , A, B are 2x2 ant i- Hermit ian trace- free matrices depending 
only on (u, v, x) and u, v, x) is an unimodular 2x2 matrix satisfying the 
reality condition: %l)(\)ip(X)^ = I. The gauge choice of the matrices A, B can 
be obtained by setting for A = 0: ^(A = 0) = J(u, v, x)~ x . Then the system 
(J2EJ) implies that 

A = J- l d u J, B = r x d x J. (27) 

Using the standard method of Riemann-Hilbert problem with zeros, Ward 
jU] was able to construct the multi-soliton solutions of (J23J) by assuming that 
ijj has a simple poles in A and no others, ie: 

k=1 A - [L k 

where M k are 2x2 matrices independent of the spectral parameter A given 
by (|13H14j) . n is the soliton number and the complex parameter [i k defines 
the velocity of the k-th soliton. Note that the construction is similar to the 
construction of baby monopoles presented in section 3.1. Again, here the 
components of M k are given in terms of a rational function f k of the complex 
variable Uk = x + [i k u + [i^v. In particular m k a are holomorphic function 
of Uk defined as: m k a = (m^rrij) = (1,/fc). By letting to be exponential 
functions of uj k the solutions correspond to extended waves which suffer a 
phase shift upon scattering studied in In general though, by letting f k 
to be any rational function of uj k the corresponding soliton dynamics is trivial 
with no change of direction or phase shift. 

Remark: In [25 j it was shown that when the spectral theory was applied 
for (J23)) the eigenfunction satisfies a local Riemann-Hilbert problem given by: 

ijj- = ip + (i - S(x + \u + A" 1 *;)) . (29) 



and the pure soliton solutions can be obtained by letting A = pi in (jzajl while 
the components of the M k matrix can be obtained by solving simple algebraic 
matrix equations. 

4.1 Non-trivial Soliton Scattering 

Here we present the way to obtain solitons with non-trivial scattering prop- 
erties, similar to the anti-de Sitter ones. Note that since the space is flat the 
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non-trivial scattering can be observed and studied in a more clear way. Once 
more, we concentrate in cases where the formula (|28j) is not well-defined: ie 
when ji). are not distinct or when /2fc = \i\ for any (k, I). 

<0 n/N Scattering 

Assume that ip has a double pole in A and no others: 

(30) 

where are 2x2 matrices independent of A. This hypothesis can be 
generalized by taking ip to have a pole of order n in A (for more details see 
13J). Then ip satisfies the reality condition if and only if is factorized as ^2] 

where qi for i = 1 , 2 are two-dimensional vectors of the form: 

gi = (i + 1/| 2 ) (i,f) + e (//-//) (/"-l) 

Q2 = (l,f). (32) 

Here / is a rational function of u = x + fiu + n~ l v and 9 = (u — fi~ 2 v) f'(u) + 
h(u>). So the solution corresponds to a large family since one may take / and 
h to be any rational meromorphic functions of uj. Then the chiral field defined 
as: J = = 0) _1 which takes the product form 




J=ll+ v 1+ ^ ^^^M (33) 



is smooth on M 2+1 (since the vectors qi are nowhere zero) and satisfies the 
boundary condition irrespectively of the choice of / and h. 

The solution (J31|) has been obtained as a limit of the simple-pole case 
(|50|) with n = 2. The idea is to take the limit — > \i by setting \i\ = fi + e, 
fi2 = t 1 ~~ £ which implies that the corresponding functions take the form 
fi(uJi) = +eh(ui) and /2(^2) = f{vz) —eh{u)2). Note that in the limit 

e — > then ^(^2) — /i(^i) ~~ ^ which results the smoothness of r/>. The 
corresponding solitons are located when the chiral field J departs from its 
asymptotic value Jq which occurs when 9 — > 0. 
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In principal, it is possible to visualize the emerging soliton structures 
when the functions / and h are rational functions of degree p, q G IN" of 
the form: f(w) = oj p and h(u) = u q . In fact, for q > p the configuration 
consists of p — 1 static solitons located at the center-of-mass of the system 
forming a ring-like structure if more than one; accompanied bjN = q — p+l 
solitons which (initially) accelerate towards the one in the middle, scatter 
at an angle of n/N and (finally) decelerate as they separate. This follows 
from the following analytic argument: the field J departs from its asymptotic 
value Jo when 9 = lo p ^ 1 (p [u — n~ 2 v) + oo N ^j — > which holds when either 
ujP- 1 = or u N + p(u — n~ 2 v) = 0. These are approximately the locations 
of the solitons in the xy-pl&ne. 

Let us illustrate the above family of soliton solutions by examining two 
simple cases by giving specific values to the arbitrary parameters /i, f(u) and 
h{uj). 

• Let us investigate the simple case where \i = i, f{uj) = to and h(u) = 0. 
Then the corresponding energy density takes the simple form 

4 + 2 2 + 4t 2/ 2 2 + 1) + 1 

S = 16— J 9 34 

(r 4 + 2r 2 + At 2 + l) 2 K J 



where r is the polar coordinate: r = \Jx 2 + y 2 while the energy density 
is time reversible. Figure El presents snapshots of the stationary soliton 
configuration which forms a single peak at t = and expands to a ring 
as t increases. For large positive t, the height of the ring (ie maximum 
of £) is proportional to t~ l while its radius is proportional to \fi. 

Next let us take \i = i, f(w) = oj and h(u) = uj. The energy density is 
given by 

_ ^5r A + l0r 2 + 4t 2 (l + 2r 2 )-8t(x 2 -y 2 ) + l 

t = 16 (6b) 

(5r 4 + 2r 2 + At 2 + 8t(x 2 - y 2 ) + l) 2 

and is symmetric under the interchange t — > —t, x y,y x. Figure 
|7| presents the collision of two solitons near t = 0: the solitons acceler- 
ate towards each other, scatter at right angles and decelerate as they 
separate. In particular, the collision is time symmetric and elastic (no 
emitted radiation). 

Note that for negative t, the two solitons (following the asymptotic ar- 
gument of J) are located on the x-axis at x ~ ±y/t while for positive 
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Figure 6: A stationary soliton configuration at different times in flat space- 
time. 
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t, they are located in the y-axis at y « ±^/— t- Once more the corre- 
sponding solitons are not of constant size: their height is proportional 
to t^ 1 while their radii are proportional to \ft. In fact the solitons 
spread out as they move apart. 

Although it seems strange that by taking the limit of soliton solutions 
with trivial scattering soliton solutions with non-trivial scattering have been 
obtained there is an explanation as shown in [13J. By studying the effect of 
the variation of e to a two soliton configuration it has been observed that at 
the limit e —>■ the solitons disperse, shift and interact with each other since 
their internal degrees of freedom and their impact parameter change. As a 
result the two initial well-separated solitons form a ring at the limit e — > 0. 

In general in a head-on collision of N indistinguishable solitons the scat- 
tering angle of the emerging structures relative to the incoming ones can be 
n/N . In particular when the N solitons are very close together they merge 
and form a ring-like structure and finally they emerge from the ring in a 
direction that bisects the angle formed by the incoming ones. 



<0 Elastic soliton- antisoliton non-trivial scattering 



Next we construct a large family of solutions which represent soliton- 
antisoliton field configurations |13j . The corresponding solution has the form 

4 A = / + — + (36) 

A — i A + i 

where n k ,m k (for k = 1,2) are complex-valued two-dimensional vectors in- 
dependent of A of the form: 

.1 /i f\ ™2 



m 1 = (1,/), m l = (-/,!) 



n x \ 1 /2i(l + |/| 2 ) 2w \ /m lf 

n 2 ) :: (i + |/|2)2 + | w | 2 ^ -2w -2t(l + \f\ 2 )) U 2t 



(37) 



with 



w 



hf + 2tf. (38) 



So we generate a solution J = ip(X = 0) which depends on the two arbitrary 
rational functions / = f(z) and h = h(z) where z = x + iy. In the general 
case where / = z p ,h = z q for p, q > the energy of the configuration 
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is equal to: E = (2p + q) 8n and the corresponding solutions consists of 
2p+q lumps at the xy-plane which scatter non-trivially and are combinations 
of solitons and antisolitons. To prove that the corresponding configuration 
consists of solitons and antisolitons a topological charge was introduced for 
the integrable chiral model by exploiting its connections with the 0(3) sigma 
model. 

A topological charge may be defined for the chiral model (although it is 
not a topological model) by exploiting its connection with the 0(4) sigma 
model [321 133]: ie bv letting 

J = J0o + icr <fi (39) 

where cr are the usual Pauli matrices and (0o> </>) = (0o, 0i, 02, 03) is a four 
vector of real fields that are constrained to lie on S 3 with the constraint 
0o + 00 — 1- The only static finite energy solutions of the 0(4) sigma 
model correspond to the embedding of the 0(3) sigma model [31]. Therefore 
the only static finite energy solutions of ([2*3]) are the 0(3) embeddings that we 
shall describe. This is because for the one-soliton solution (static or Lorentz 
boosted in the |/-axis) the system behaves like the 0(4) model for which the 
0(3) embedding is totally geodesic. However for time-dependent configura- 
tions the evolution of the field will not general lie in an 0(3) subspace of 
0(4). 

In studying soliton-like solutions we require that the field configuration 
has finite energy. This implies that the field must take the same value at all 
points of spatial infinity so that the space is compactified from M 2 to S 2 . At 
fixed time the field is a map from S 2 into the target space. Now for the 0(3) 
model the field is a map : S 2 — > S 2 and due to the homotopy relation 

n 2 (S 2 ) = Z (40) 

such maps are classified by an integer winding number Af which is a conserved 
topological charge. An expression for this charge is given by 

N = -!- ( eij 0> {di<t> A djtj>) d 2 x (41) 

S7T J 

where % = 1, 2 while x l = (x, y). 

On the contrary, for the 0(4) sigma model the field at fixed time is a map 
(0o, 0) : S 2 — > S 3 which implies that 

7r 2 (S 3 ) = (42) 
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and therefore there is no winding number. However, for soliton solutions 
that correspond to some initial embedding of 0(3) space into 0(4) there is 
a useful quantity presented below. 

Consider the 0(4) configuration which at some time corresponds to an 
0(3) embedding which we choose to be <po = for definiteness. At this 
time the field is restricted to an S 2 equator of the possible S 3 target space. 
Suppose that the field never maps to the antipodal points (1,-1) so the 
target space is Sq ~ S 2 x R and thus we have the homotopy relation 

vr 2 (5 3 ) = tt 2 (S 2 x R) = ti 2 (S 2 ) © vr 2 (R) = Z (43) 

and therefore a topological winding number exists. An expression for this 
winding number is easy to give since it is the winding number of the map 
after projection onto the chosen S 2 equator ie 

where 4>' = <p/\<p\. If the field does map to the antipodal points at some time 
the winding number is ill defined at this time and if considered as a function 
of time M' will be integer valued but may suffer discontinuous jumps as the 
field moves through the antipodal points. In the following examples before 
comparing the solution J with the 0(3) embedding it is convenient to perform 
the transformation J — > MJ with 

M ^7i(-i \) (45) 

so that the evolution of the field remains close to the 0(3) embedding. 
One particular example is given below: 

• Lets take the simplest case where f(z) = z and h(z) = 1. Roughly 
speaking the chiral field J departs from its asymptotic value J when 
w — > which implies that z — > ±\/—2t. More precisely, the two 
structures, for t negative, are located approximately at x ~ ±a/— 2t 
while, for t positive, they are on the y-axis at y ~ ±v2t. Figure |H1 
illustrates the scattering behavior of our configuration near t = 0. 

The picture is consistent with the properties of the energy density of 
the solution which is given by 

2r 4 + 4r 2 + 4t 2 (l + 2r 2 )-4t(x 2 -y 2 ) + l . . 

£ = 16 y — i y \ 46 

(2r 4 + 2r 2 + At{x 2 - y 2 ) + At 2 + l) 2 
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Figure 8: An elastic non-trivial scattering of a soliton-antisoliton configura- 
tion in flat spacetime. 
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Figure 9: Topological charge density at increasing times for soliton- 
antisoliton elastic scattering. 



and is symmetric under the interchange t — > — t, x — > y, y — > x. The 
time symmetry of the energy density confirms the lack of radiation. 
However, the corresponding localized structures are not constant of 
size; in fact, their height is proportional to t^ 1 and their radius is 
proportional to \fi. 

The projected topological charge H' is zero throughout the scattering 
process while the topological density has an almost identical distribu- 
tion (up to a scale) to that of the energy density as shown in Figure 
HJ Therefore, the configuration represents a soliton and an antisoli- 
ton that are clearly visible as distinct structures having respectively 
one and minus one units of topological charge concentrated in a single 
lump. 

In the general case where h = z p and h = z q , the chiral field J departs 
from its asymptotic value Jo when w = z (j> ~ 1 \2tp+z N ) — ► with N = p+q+1 
which is true when either z^" 1 ' = or 2tp+z N = and this is approximately 
where the lumps are located. So J represent a family of soliton-antisoliton 
solution which consists of (p — 1) static soliton-like objects at the origin with 
iV others accelerating toward them, scattering at an angle of tt/N and then 
decelerating as they separate. 

Note that in the non-integrable models (usually) there is an attractive 
force between solitons of opposite topological charge which implies that the 
initially well separated solitons and antisolitons attract each other and even- 
tually annihilate into a wave of pure radiation which spreads with the velocity 
of light in a direction perpendicular to the motion of the initial structures 
as shown in [2E]- However, it is known that the interaction forces between 
solitons and antisolitons do depend on their relative orientation in the in- 
ternal space which implies that the cross section for the soliton-antisoliton 
elastic scattering is different than zero. In fact the proton-antiproton elastic 
scattering is seen in a reasonable fraction of cases. Therefore the analytic 
construction of families of soliton-antisoliton configurations with elastic scat- 
tering properties obtained are of immense importance since not only provide 
a major link between integrable and non-integrable models but can actually 
have physical applications in the laboratories. 
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5 Conclusions 



The infinite number of conservation laws associated with a given system place 
severe constraints upon possible soliton dynamics. The construction of exact 
analytic multi-solitons with trivial scattering properties is a result of such 
integrability properties. In this paper soliton and soliton-antisoliton struc- 
tures have been presented for two planar integrable equations related with 
the three-spatial dimensional self-dual Yang-Mills-Higgs equations. These 
structures travel with non-constant velocity, their size is non-constant and 
they interact non-trivially. As we have already mentioned this non-trivial 
scattering is not usual in an integrable theory but is exceptional. Such 
results are useful for investigating the connection between integrable and 
non-integrable systems which possess multi-soliton solutions. In addition, 
they indicate the likely occurrence of new phenomena in higher-dimensional 
soliton theory which are not present in one- dimensional systems. 

The multi-pole ansatz have been extended to other planar integrable sys- 
tems in order to obtain soliton solutions with non-trivial dynamics. In partic- 
ular, a large class of solutions of the Davey-Stewartson II equation has been 
constructed in [33] which have an arbitrary rational localization in the plane 
and describe typical interactions consisting of head-on collisions with a 71/ 2 
scattering angle. Similarly, in [36| the discrete spectrum of the non-stationary 
Schrodinger equation and localized solutions of the Kadomtsev-Petviashvili 
I equation are studied via the inverse scattering transform. It was shown 
that there exist infinitely many real and rationally decaying potentials which 
correspond to a discrete spectrum whose related eigenfunctions have multi- 
ple poles in the spectral parameter; while the resulting localized solutions 
of the Kadomtsev-Petviashvili I behave as a collection of individual humps 
with nonuniform dynamics. 

It would be interesting to understand the role of higher poles in algebraic- 
geometry approach like twistor theory (for example, the function \I/ given by 
f!16|) correspond to n = 2 bundles), and also to investigate the construction of 
the corresponding solutions and their dynamics in de Sitter space. Finally, it 
would be interesting to extend our construction in higher dimensional gauged 
theories and investigate the scattering behavior of the corresponding classical 
solutions and, also, consider and study its noncommutative version (see, for 
example, Ref. [37]). 
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